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ABSTRACT 

A  long-period  forecast  is  made  utilizing  a  two-level  quasi- 
geostrophic  model.   The  model  includes  friction  and  heating  which  is 
a  linear  function  of  y.   The  model  was  further  simplified  by  restrict- 
ing the  disturbance  to  one  wave  in  the  zonal  direction.   Experiments 
were  performed  with  two  distances  between  the  walls.   In  the  case  of 
the  longer  separation,  a  solution  with  appreciable  time  fluctuations 
was  obtained  with  the  largest  fluctuation  being  investigated  in  more 
detail.   The  smaller  separation  revealed  a  traveling  wave  of  constant 
amplitude  similar  to  that  observed  in  certain  dishpan  experiments. 
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I.   INTRODUCTION 

Knowledge  of  the  general  circulation  has  been  made  possible  with 
the  advent  of  high-speed  large-capacity  electronic  digital  computers. 
In  a  classic  experiment,  Phillips  [1956]  performed  a  long-period 
numerical  forecast  with  a  two-level  quasi-geostrophic  model  in  a 
zonally  periodic  channel  on  a  3  plane.   Friction  and  heating  were 
included  in  his  equations  with  the  latter  as  a  linear  function  of 
latitude.   He  introduced  a  finite  amplitude  initial  disturbance  into 
a  baroclinically  unstable  zonal  current.   Major  deficiencies  were  the 
error  from  the  3  plane  and  the  geostrophic  approximations,  the  latter 
making  it  impossible  to  account  for  non-geostrophic  dynamics.   Further- 
more, due  to  computational  errors,  he  did  not  achieve  a  statistically 
steady  state.   However,  Phillips  did  initially  demonstrate  the  feasi- 
bility of  numerical  simulation  of  the  atmosphere. 

General  circulation  experiments  have  progressed  in  sophistication 
and  complexity  since  Phillips'.   Smagorinsky  [1963]  formulated  a  two- 
level  baroclinic  primitive  equation  model  by  using  the  Eulerian 
equations  of  motion  within  a  spherically  zonal  strip.   Continued 
sophistication  of  the  experiments  resulted  in  a  nine-level  model 
utilizing  the  primitive  equations  of  motion  [Smagorinsky,  Manabe 
and  Holloway,  1965],   This  study  included  a  resolution  of  boundary 
layer  fluxes  and  radiative  transfer  involving  ozone,  carbon  dioxide 
and  water  vapor.   A  continuation  of  this  study  was  performed  for  the 
tropics  which  was  successful  in  simulating  the  tropical  convergence 
zone  and  the  thermal  structure  [Manabe  and  Smagorinsky,  1967].   Mintz 


[1964]  conducted  an  experiment  in  which  the  primitive  equations  of  motion 
were  developed  for  a  two-level  model  of  the  atmosphere  with  the  heating 
and  friction  terms  retained.   This  model  further  brought  in  the  effects 
of  land-sea  contrasts  in  heating.   Nitta  [1967]  utilized  a  twenty-level 
model  for  computing  vertical  distribution  of  the  geopotential  flux  and 
conversion  of  eddy  potential  energy  and  eddy  kinetic  energy  due  to  non- 
stationary  disturbances.   A  model  which  has  the  distinguishing  feature 
of  using  height  rather  than  pressure  was  developed  [Kasahara  and 
Washington,  1967]  and  hydrostatic  equilibrium  was  maintained  in  the 
system.   Some  advantages  of  this  model  are  the  prognostic  equations 
have  a  simpler  form  than  those  of  the  a  system  and  the  lower  boundary 
conditions  may  be  easily  formulated.   The  integration  of  the  nine-level 
primitive  equation  model  was  further  extended  by  the  box  method  [Kurihara 
and  Holloway,  1967]  which  conserved  the  conservational  properties  of  the 
original  equations.   Recently,  the  effect  of  hydrology  of  the  earth's 
surface  was  incorporated  in  the  nine-level  primitive  equation  model 
[Manabe ,  1969].   The  scheme  involved  the  prediction  of  water  vapor  in 
the  atmosphere  and  the  prediction  of  soil  moisture  and  snow  cover.   The 
numerical  integrations  were  performed  for  the  annual  mean  distribution 
of  solar  insolation. 

The  increased  sophistication  and  complexity  of  recent  models  have 
also  produced  complex  results  which  make  interpretation  difficult. 
Therefore,  in  this  study,  the  complexities  have  been  kept  to  a  minimum 
while  still  maintaining  a  description  of  the  basic  dynamic  processes. 
When  the  effects  of  heating  and  friction  are  included,  the  two-level 
quas i-geostrophic  model  is  the  simplest  model  capable  of  describing  the 
general  circulation. 


The  model  is  further  simplified  in  that  the  disturbance  is  restricted 
to  a  single  wave  in  the  zonal  direction.   Phillips'  [1956]  experiments 
and  Smagorinsky 's  [1963]  experiment  which  contained  no  forcing  function 
in  the  east-west  direction,  displayed  a  single  predominate  wave  number 
which  contained  much  of  the  disturbance  energy.   The  y-structure  of  the 
disturbance  and  the  mean  flow  are  calculated  accurately  with  a  sufficient 
number  of  grid  points.   Some  general  circulation  experiments  have  been 
made  with  even  simpler  models  which  represent  the  y-variation  of  all 
quantities  with  a  few  terms  of  a  Fourier  series  [Bryan,  1959;  Lorenz, 
1962,  1963;  Young,  1966],   These  models  do  not  provide  an  adequate 
description  of  the  barotropic  interaction  between  the  disturbance  and 
the  mean  flow.   Consequently,  this  model  appears  to  be  the  simplest 
model  which  can  describe  the  main  mechanisms  of  the  general  circulation 
of  the  atmosphere. 


II.      THE   FORECAST   EQUATIONS 


A  simple    two-level  model    is   constructed  by  dividing   the   entire 
atmosphere    into   four    layers    of  constant    pressure   differential,  -r*- 
(fig.    1),    numbered   0    to  4    from   top    to  bottom.      Assume   vertical  motion 
to  be   zero  at   the   top   of  the   atmosphere,   while 
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Fig.  1.   Two-level  model  used  for  prediction 

the  vertical  motion  term  at  the  earth's  surface  is  approximated  by 

cd  =  -pgw4>  (2.1) 

Charney  and  Eliassen  [1949]  used  the  Ekman  theory  to  derive  an 
expression  for  w,  which  is 

i  r2  V1/2 

w4  =  2  CfTV   sin  2  cc  C4>  (2.2) 


where  A   is  the  vertical  eddy  viscosity  and  a,  the  inflow  angle.   The 
surface  geostrophic  vorticity  is  approximated  by 

£4=C3=V2T3.  (2.3) 

Following  the  development  of  Thompson  [1961],  but  with  different 
notation,  begin  with  the  quasi-geostrophic  vorticity  equation 

|-  A  +  k   x  VyV(A)  +  eo  §f  -  fQ  |2  =  AH  A.    (2.4) 
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Apply    this   equation   at    levels    1   and   3   giving 


dY, 


uv 


!?  A +  lk x *vv<*2v +  po i1  -  fo i  ■  ah v\ 


(2.5) 


dY. 


-  V  Y3  +  ,k  x  VY3.V(V  Y3)  +  0Q  ^-  +  fQ-^—  =  ^  V  y3.         (2.6) 


^t 


If  it  is  desired  to  stop  at  some  intermediate  point,  say  the  tropopause , 
smaller  layers  could  be  utilized  rather  thau  to  include  the  entire  atmosphere 

Next  consider  the  quasi-geostrophic  first  law  of  thermodynamics  in 
the  form 


iL  $1 

dt  dp 


Ik  x  VY-V  (Ip  +  f-  a)  = 


V 


where 
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R  T  rdT 
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(2.7) 


(2.7a) 


and  T   is  the  temperature  from  the  standard  atmosphere.   When  this 


equation  is  applied  at  level  2,  the  result  is 


(W 


3t  <VV+lkxV"  -2 


V(Y1-Y3) 


AP  au>2       H 
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(2.8) 


Define  the  following  quantities 
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- 
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(2.9) 
(2.10) 


which    implies 


*1  =  \  +  V 

t3  =  ■■»„,»  rT. 


(2.11) 
(2.12) 


Here  YT  is  proportional  to  the  layer  thickness  and  is  therefore  a 
measure  of  the  mean  temperature.   Using  these  definitions,  add  (2.5) 
and  (2.6)  and  divide  the  result  by  2,  obtaining 
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n         2  2  2  ^  ^A 

%r  V  Y     +  !k  x  VY    -V(V  Y)  +  k  x   YT-V(V  YT)  +  0     XT  "    f     T7T 
dt  m  m  m  1  T  o    dx  o2Ap 


=  AH 


V  Y 


m 


(2.13) 


For  the  second  forecast  equation,  subtract  (2.6)  from  (2.5)  and 
eliminate  u)_  usi,ng  (2.8) 

|r  (V2V)  \  +  Ik  x  VYm-V(V2-n2)  YT  +  k  x  VYT-V(V2Ym) 


at 


+  3  -si  +  f  ~- 

o  dx    o  2Ap 


2  KQ2       4 
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(2.14) 


where 


9    2  f 
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|i  -   2 

Ap  a 


(2.14a) 


These  are  the  prediction  equations  for  the  model.  At  this  point 
we  restrict  ourselves  to  a  disturbance  of  wave  number  k  and  the  mean 
flow.   The  fields  may  be  defined  as  follows 

Y  =  E(y,t)  +  A(y,t)  cos  kx  +  B(y,t)  sin  kx,     (2.15) 

m 

YT   =  F(y,t)  +  C(y,t)    cos   kx  +  D(y,t)    sin  kx ,  (2.16) 

where   k   is    the  x  wave   number,   A  through  D  are   Fourier  amplitudes    of   the 
disturbance   and   E  and  F  are   the   zonal  mean   fields.      Now  substitute 


expressions    for   Y     and   YT   into    (2.13),    separate    the  various    terms, 
neglecting   all    terms   with  wave    number   2k   or   higher.      Equating   coeffici- 
ents   of   the   cosine   terms   gives 


3_fjTA 
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For   the   sine   terms    the   result    is 


(H  -  «0  - 
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(2.18) 


Repeating   the    procedures    for    (2.14),    the   equations   are 
2        „   2 


d     f''b  C 
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and 
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dy  dy 


(2.20) 
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Equating  coefficients  of  terms  independent  of  x  in  (2.13)  gives 


a   a  e  =  k  a  r    a'B 
"  2  ay  L    a  2 


B   *LA  +  C  4   "   D   ^ 
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3, 


ay 


2J 
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(2.21) 


Similar ly,  coefficients  indepdendent  of  x  in  (2.14)  results  in 

2^  M.  _  is  a_ 
y  ay     2  ay 


^y       2^  ap     k  a_ 

^2  "  ^ay  =  2  ay 


ay        ay        cv        ay 
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K  (M  .  MT)  .  ^    K  Q2  +  A    cl  firs 

K  vay    *yJ       2       +  ah  ay  ^3> 


(2.22) 


where 


K  = 


fg/A 
o     /    m 

2RT      fd 


1/2 


(2.23) 

The  prediction  equations  are  now  in  the  form  to  be  used  in  the  model. 
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III.   FINITE  DIFFERENCE  SCHEME  AND  BOUNDARY  CONDITIONS 

The  finite  difference  scheme  used  is  illustrated  below  with  a 
sample  variable  N; 

2 
'  -  K  (Y.  .,  -  2N.  +  N ._/)  (3.2) 


Sy2   H2Vi+1 

dy    2H 


"  2Ni+l  +  O  -  (Ni  -  2Ni-l  +  Ni-2)]    (3'3> 


where  i  is  the  grid  index  and  H  the  distance  between  grid  points. 

Centered  time  differences  are  used  for  all  quantities  except  those 
involving  friction.   The  frictional  terms  are  computed  at  time  (t-At). 
In  all  cases  the  first  step  was  a  forward  time  step.   A  special  time 
step  was  introduced  every  24  time  steps  to  avoid  the  separation  of 
solutions  at  even  and  odd  time  steps.   These  special  time  steps  used 
the  backwards  difference  method  utilized  by  Matsuno  [1966].   Second 
order  equations  (2.17),  (2.18),  (2.19)  and  (2.20)  for  time  tendencies 
are  solved  by  the  exact  method  of  Richtmyer  [1957].   A  modification  of 
the  Richtmyer  method  used  by  Phillips  [1956]  was  necessary  for  equation 
(2.22)  to  satisfy  the  boundary  conditions.   Equation  (2.21)  was  solved 
by  a  direct  marching  process  utilizing  the  boundary  conditions. 

The  following  boundary  conditions  were  established: 

The  quantities  A  through  D  =  0 

d  dF    _    .  d  dE  , ,        _ 

■v—  -n-  =  0  and  v~  -T-  -0      all  at  y  =  0  and  y  =  W. 


15 


Following  Phillips  [1956]  the  mean  and  disturbance  vorticity  are 
set  equal  to  zero  at  the  boundary  in  order  to  make  the  boundary  walls 
smooth. 
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IV.   INITIAL  CONDITIONS 

A  natural  initial  state  would  be  one  of  constant  temperature  and 
no  motion.   The  heating,  in  time,  would  build  up  a  baroclinic  current. 
No  disturbance  could  develop  until  the  vertical  shear  of  the  zonal 
current  reached  a  critical  value.   In  order  to  save  computer  time,  the 
initial  state  was  chosen  to  consist  of  a  zonal  current  with  a  vertical 
shear  near  the  critical  value  [see  Thompson,  1961].   The  initial  dis- 
turbance has  a  wave  length  and  a  y-structure  such  that  when  superimposed 
on  a  baroclinic  current  with  no  horizontal  shear  it  would  have  maximum 
growth  rate . 

The  initial  conditions  for  the  model  are 


A  =  200  sin  7^ 
w 


B  through  D  =  0 

f  HW6 

r-    o  ny 

F  =  —  cos  rf- 

TT  W 

f  HW6 

E  = cos  -f- 

TT  W 
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V.   RESULTS 

The  main  experiments  were  conducted  using  the  physical  parameters 
listed  in  Table  I. 


f 

o 

= 

-4 
1  x  10   per  second 

a 

= 

2.3  meters  squared  per  second  squared  per  centibar 
squared 

0 

= 

1.67  x  10    per  meter  per  second 

AP 

= 

50  centibars 

A 

V 

= 

40  meters  squared  per  second 

AH 

= 

1  x  10  meters  squared  per  second 

a 

= 

22.5  degrees 

2 

= 

-12 
2.49  x  10    per  meter  squared 

H 

= 

200  kilometers 

At 

= 

0.5  hours 

k 

= 

2tt/4000  per  kilometer 

Table  I  -  Numerical  Values  of  the  Physical  Constants 
The  variation  of  experiments  2  and  3  from  the  main  experiment  is 
tabulated  in  Table  II. 


W  (Km) 

Q  (Kj  per  ton  per  sec) 

Experiment  1 
Experiment  2 
Experiment  3 

8000 
8000 
4000 

4  x  10"3  (y  -  £)/w 
8  x  10'3  (y    -   D/W 
8  x  10"3  (y   -  D/W 

Table  II  -  Variation  of  the  Physical  Parameters 
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The  forecast  period  for  all  three  experiments  was  300  days.   Figure 
2  depicts  a  plot  of  the  disturbance  energy  for  the  300  day  forecast 
period.   The  disturbance  energy  referred  to  throughout  the  discussion 
is  defined  in  Appendix  C  (Eqs .  C-10  and  C-ll).   A  main  feature  of  Fig. 
2  is  that  the  disturbance  energy  does  not  reach  a  steady  state  although 
it  appears  to  be  statistically  stable  after  the  first  120  days. 

Figure  3  displays  a  time  section  of  the  mean  wind  for  the  period 
120-260  days  which  is  the  main  time  period  of  interest  for  this  experi- 
ment.  The  general  increase  of  the  disturbance  energy  begins  to  level 
off  near  day  140  which  corresponds  to  the  initial  appearance  of  the 
southern  jet.   This  result  is  similar  to  observations  of  Riehl,  Yeh 
and  LaSeur  [1950].   The  period  between  day  190  and  day  200  is  of  con- 
siderable interest.   The  disturbance  energy  has  its  greatest  fluctuations 
during  this  period  and  it  corresponds  to  the  termination  of  the  southern 
jet. 

Figures  4-7  represent  the  amplitudes  and  phases  of  Y  ,  V  and  w~ 
in  the  form  Pcos(kx-S)  where  P  is  the  amplitude  and  6  is  the  phase. 
Figures  4  and  5  represent  day  190,  prior  to  the  disturbance  peak,  and 
Figs.  6  and  7  correspond  to  day  196,  after  the  peak.   The  necessary 
condition  for  barotropic  instability  is  that  g  -  — — tt~   change  sign  some- 
where  in  the  region  [Kuo,  1949].   Calculation  at  day  190  showed  that 
this  condition  was  satisfied  in  the  vicinity  of  the  southern  jet.   It 
is  seen  in  the  lower  portion  of  Fig.  4,  between  y  =  1000  km  and  y  = 
2500  km,  that  the  phase  tilts  of  Y  ,  Y„  and  w2  are  opposite  to  the 
shear  of  u..  .   These  phase  relationships  indicate  barotropic  amplifi- 
cation and  therefore  a  transformation  of  energy  from  the  mean  flow  to 
the  disturbance.   However,  the  northern  portion  of  Fig.  4  shows  the 
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same  tilt  as  the  shear  of  u  which  indicates  barotropic  damping:   this 

implies  an  energy  transformation  of  disturbance  kinetic  energy  to  mean 

kinetic  energy.   Figure  5  portrays  the  amplitude  of  Y  and  Y_ ;  clearly 

m      T 

Y  is  larger  and  both  fields  have  a  mid- latitude  maximum.   A  secondary 

peak  of  amplitude  in  Y   is  beginning  to  form  which  is  a  product  of  the 

barotropic  instability.   Figure  6  shows  that  all  parameters  O  .  Y  .  w  ) 

v  v  m'   T"   2y 

have  a  phase  tilt  which  indicates  barotropic  damping  over  the  entire 
range.   The  amplitude  of  Y  and  Y  corresponding  to  Fig.  6  are  displayed 
in  Fig.  7.   Of  primary  note  is  the  distinct  double  maximum  of  amplitude 
which  is  a  result  of  barotropic  instability  in  the  southern  region. 
Figure  8  depicts  the  energy  transformations  over  a  period  which  includes 
the  maximum  disturbance  energy.   The  first  item  of  note  is  the  nearly 
exact  correlation  between  [P-P']  and  [P1 «K' ]  which  in  an  energy  balance 
would  almost  completely  offset  one  another  thereby  leaving  a  small 
amount  of  disturbance  potential  energy  change.   It  should  be  noted 
that  just  before  the  event  of  maximum  disturbance  energy,  the  barotropic 
damping  [K' -K]  approached  a  value  of  zero.   This  fact  permits  the  dis- 
turbance kinetic  energy  to  increase  at  a  rapid  rate  due  to  the  large 
value  of  [P'-K1]  and  to  continue  to  increase  until  [k'-K]  attains  a 
value  large  enough  to  terminate  growth. 

Figure  9  depicts  the  disturbance  energy  over  the  forecast  period 
for  experiment  2.   As  stated  earlier  (Table  II),  this  experiment  has 
twice  the  heating  rate  of  experiment  1.   A  steady  state  was  not 
attained  but  the  disturbance  energy  became  statistically  steady.   The 
main  differences  from  experiment  1  are  that  the  disturbance  energy 
fluctuates  more  rapidly  and  that  its  magnitude  is  much  higher. 
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A  depiction  of  the  disturbance  energy  over  the  forecast  period  for 

experiment  3  is  contained  in  Fig.  10.   This  experiment  differs  from 

experiment  1  in  that  the  disturbance  between  the  walls  has  been  cut  in 

half.   The  heating  coefficient  was  doubled  from  its  value  in  experiment 

1.   The  disturbance  energy  reached  a  steady  state  condition  shortly 

after  day  120  and  stayed  relatively  constant  to  the  end  of  the  period. 

Thus  we  have  a  single  wave  traveling  through  the  atmosphere  with  no 

variation  in  amplitude.   This  result  is  similar  to  certain  dishpan 

experiments  with  an  inner  core  [Hide,  1953;  Fultz,  et  al,  1959]. 

Figure  11  represents  the  phases  of  Y  ,  L  and  w  at  day  200  and  reveals 

that  the  parameters  have  similar  phase  variations.   The  phase  tilts 

indicate  barotropic  damping  since  there  is  a  single  jet  in  u1  (see 

Fig.  13).   Since  Y™  lags  Y  ,  there  is  a  northward  heat  flux  and  energy 

1       m 

transformation  from  the  mean  potential  energy  to  the  disturbance 

potential  energy.   The  correlation  mentioned  above  gives  a  transformation 

of  disturbance  potential  energy  to  disturbance  kinetic  energy  since  there 

is  a  tendency  for  cold  air  to  sink  and  warm  air  to  rise.   Figure  12 

represents  a  plot  of  the  amplitudes  of  Y  and  Y~  at  day  200.   The 

amplitude  of  Y   is  greater  than  Y™  and  both  have  a  single  maximum, 
m  T 

Figure  13  depicts  the  relationship  between  u„ ,  which  could  be  con- 
sidered the  surface  wind,  and  u1  at  day  200.   This  mean  wind  field  is 
similar  to  that  observed  in  the  atmosphere,  although  the  gradients  are 
much  less.   At  the  surface,  u  is  westerly  in  the  mid  latitudes  and 
easterly  at  high  and  low  latitudes  while  at  the  upper  levels  the  flow 
is  westerly  with  a  mid-latitude  jet.   Figure  14  is  a  plot  of  the 
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amplitude  of  the  mean  meridional  component  of  the  wind  (v,)  for  day  200. 
The  amplitude  of  v  reveals  a  northerly  component  at  mid  latitudes  and 
a  southerly  component  at  high  and  low  latitudes  which  implies  an  indirect 
Ferrel  circulation  at  mid  latitudes  and  a  direct  Hadley  circulation  at 
high  and  low  latitudes, 
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VI.   CONCLUSIONS 

A  two-level  quasi-geostrophic  model  was  formulated  for  simplified 
studies  of  the  general  circulation  of  the  atmosphere.   The  model  was 
further  simplified  by  restricting  the  disturbance  to  one  wave  in  the 
zonal  direction. 

Numerical  integration  was  computationally  stable  for  a  period  of  at 
least  300  days.   In  all  the  experiments  the  same  heating  function  was 
used  which  was  linear  in  y  and  independent  of  time.   The  main  experi- 
ment with  a  wall  separation  of  8000  km  produced  fluctuating  solutions 
with  irregular  time  variations.   Attention  was  focused  on  a  sub  period 
of  the  integration  where  the  disturbance  featured  a  rapid  growth  followed 
by  a  rapid  decay.   Apparently  a  double  jet  structure  in  the  mean  flow 
became  barotropically  Unstable  at  the  beginning  of  the  period  which 
resulted  in  the  rapid  growth  of  the  disturbance.   By  the  end  of  the 
period,  a  single  jet  had  replaced  the  double  jet  structure.   Energy 
transformations  verify  the  roll  of  [K1 -K]  .in  the  evolution  of  this 
disturbance . 

Experiment  2  featured  a  double  heating  rate  and  results  were 
similar  to  experiment  1  with  increased  fluctuations  and  higher  energy 
values  of  the  disturbance  energy. 

Experiment  3  differed  from  experiment  1  in  that  the  distance 
between  walls  was  reduced  to  4000  km  and  the  heating  rate  was  doubled. 
The  disturbance  energy  became  constant  after  approximately  120  days 
and  remained  so  throughout  the  forecast  period.   This  situation 
represents  a  wave  of  constant  amplitude  propagating  through  the 
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atmosphere.   The  behavior  resembles  the  flow  observed  in  certain  dish- 
pan  experiments  which  contained  a  central  core  [Hide,  1953;  Fultz ,  et  al, 
1959].   Apparently  the  explanation  of  the  difference  in  behavior  between 
experiment  1  and  experiment  3  is  related  to  the  fact  that  for  W  =  8000  km. 
at  least  two  energy  producing  modes  of  disturbance  energy  can  exist  in  y, 
whereas  for  W  =  4000  km  only  one  energy  producing  mode  can  exist.   The 
presence  of  two  energy  producing  modes  apparently  leads  to  fluctuations 
in  time  but  the  one  modal  case  produces  a  steady  wave. 

Further  studies  utilizing  this  simple  model  could  be  performed  such 
as  longer  period  forecasting;  changes  in  the  form  of  heating;  simulating 
dishpan  experiments  by  setting  3  "  0  and  concentrating  the  heating  near 
the  boundaries,  and  a  change  in  the  wave  number  k  to  name  a  few.   Any 
or  all  of  these  studies  could  further  enhance  our  understanding  of  the 
general  circulation  of  the  atmosphere, 


37 


APPENDIX  A  -  COMPUTATION  OF  VERTICAL  MOTION 


The  computation  of  vertical  motion  (ax, )  is  necessary  for  the  energy 
equation.   To  derive  the  computational  equation  for  vertical  motion, 
start  with 


r^Yn 


^2   ^  Ap^ 


dl 


^  HQ2- 

+     U         -ST +     V         -=7 +    —T— 

_ot  m   ox  m   dy  I    . 


(A-l) 


where 


i  ^ 

1         m 

m            f  ay 

o  J 

1       m 

m        f      ax 
o 


(A-2) 


(A-3) 


Substituting  equations  (A-2),  (A-3),  (2.15)  and  (2.16)  into  (A-l), 
collecting  terms,  and  again  neglecting  terms  with  wave  number  2k  or 
higher  results  in 


OF 


2       J         .      ToC       k     /oE  „ 
^2  =  Wo  f  °S  **  Ldl  "  T  ^  D  "  o^  ] 


sin  kx 


_dB        k_  rdE  dF 

^t  +  f    Voy  C   "  oy 


0J*£ 


,       *Q2( 


(A-4) 


Both  the  disturbance  and  the  average  vertical  motion  are  given  by  this 
equation,  which  is  the  equation  used  in  the  model  for  the  computation 
of  vertical  motion. 
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APPENDIX  E  -  COMPUTATION  OF  AVERAGE  NORTH-SOUTH  WIND 


The  computation  of  the  average  north-south  wind  (v..)  is  necessary 
for  the  energy  calculations  .   The  computation  for  v..  can  be  derived 


from  Fig.  B-l  by  the  following 

+  vMV*  -   0  (B-l) 


a -y- 


CO  ———>       co   <—— —  Od 

vl       I       vl  A 


A  A   V   A   A 


CO-    ov 


2  +   1  _  0  (B.3)    ________________ _ 

Ap    dy  Fig.  B-l 

In  finite  difference  form  equation  B-3  becomes 

r*       r<*      _r("2)i+  s2)i+i>i 

(Vl   S    <Vi    "  Ip\  2  V  (B"4) 


with  v.    =  0.0 
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APPENDIX  C  -  ENERGY  TRANSFORMATION  EQUATIONS 


The  following  definitions  are  utilized  in  the  energy  transformation 
equations  of  Phillips  [1956] 


2  r.Y,  -  Y, 


Mean  Potential  Energy   P  -  ^-J  l^— — ^ J     ^ 


)isturbance  Potential  Ene         =  fej^'i  "  Y'3)"  dy 


2 

l_ 

4W 


~   1   f  -  2   —  2 
Mean  Kinetic  Energy   K  -  77:  J  (u   +  u_  ) 


dy 


(C-D 


(C-2) 


(C-3) 


^2  M,-}    ^y'  2  M*    2 


Disturbance  Kinetic  Energy   K'  -  ^{[("^t)  +  Gy^)  +  ("St)  +  C"5r)  1  dy 

(C-4) 
Phillips  [1956]  derived  the  following  disturbance  energy  trans- 
formation equations  which  provide  an  excellent  means  for  examining 
the  extent  to  which  the  model  contains  the  physical  processes  known 
to  be  important  in  the  atmosphere. 
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The  symbolic  notation  of  the  form  [P'P1]  for  example,  signifies  a 
transformation  of  energy  from  one  form  -  the  first  in  the  bracket  -  to 
the  second  form.   The  example  [p*P*]  represents  the  transformation  of 
mean  potential  energy  into  the  disturbance  potential  energy. 

Utilizing  equations  (2.9),  (2.10),  (2.15)  and  (2.16),  the  components 
of  the  disturbance  energy  become 
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Similarly,  equations  (C-5) through  (C-9)  become 
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[K'-V=lSr     /    V  [(A-C)2  +   (B-D)2]  + 
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(C-16) 
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